n=64;
w=2/3;
figure('units', 'normalized', 'position', [0, 0, 0.5, 1]);
nu = [0, 0;
      0, 2;
      1, 1;
      2, 0;
      2, 2;
      4, 0];
for j = 1 : 6
    clf;
    axis([0, 64, 0, 1])
    for k = 1:32
        lambda_1 = 1-2*w*(sin(k*pi/2/n))^2;
        lambda_2 = 1-2*w*(sin((n - k)*pi/2/n))^2;
        C1(k) = lambda_1^(nu(j, 1) + nu(j, 2))*(sin(k*pi/2/n))^2;
        C1(n-k) = lambda_1^(nu(j, 1))*lambda_2^(nu(j ,2))*(sin(k*pi/2/n))^2;
        C2(k) = lambda_2^(nu(j, 1))*lambda_1^(nu(j, 2))*(cos(k*pi/n/2))^2;
        C2(n-k) = lambda_2^(nu(j, 1)+nu(j, 2))*(cos(k*pi/2/n))^2;
    end
    plot(1:1:32, C1(1:32), '-r');
    hold on; 
    plot(32:1:63, C1(32:63), '-b');
    plot(1:1:32, C2(1:32), '-g');
    plot(32:1:63, C2(32: 63), '-c');
    lgd=legend('c1', 'c2', 'c3', 'c4');
    set(lgd, 'Fontsize', 20);
    set(gca, 'Fontsize', 20);
    s = ['../Exercise9_47_', num2str(j), '.eps'];
    print(gcf, '-depsc', s);
    hold off;
end